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Abstract 



We derive the exact critical couplings (x*,?/*), where y^jx* = y 1 + \/2 = 

1.533..., for the polymer adsorption transition on the honeycomb lattice, 

along with the universal critical exponents, from the Bethe Ansatz solution 

of the 0(n) loop model at the special transition. Our result for the thermal 

scaling dimension, and thus the crossover exponent (/> = |, is in agreement 

with an earlier result based on conformal invariance arguments. Our result 

for the geometric scaling dimensions confirms recent conjectures that they are 

given by /i^+i,3 in the Kac formula. 
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A long flexible polymer in a good solvent with an attractive short-range force between 
the polymer and the container wall is known to undergo an adsorption transition llj-Q. 
A standard model for this phenomenon is a self-avoiding walk (SAW) on a rf-dimensional 
lattice interacting with a (rf — l)-dimensional substrate. In the lattice model, the SAW has 
a Boltzmann weight x per monomer (in the bulk), with weight y per adsorbed monomer (on 
the substrate). At the adsorption transition, y*, the number of adsorbed monomers scales 
with the total length L as La ~ i^"^, where is a crossover exponent. The polymer is in the 
adsorbed phase for y > yl and the desorbed phase for y < yl, where the surface attractions 
are not effective. In the language of surface critical phenomena, the adsorption transition is 
a special transition ^ . 

Two-dimensional polymers are not without experimental interest and the above model 
has been widely studied via a number of techniques (see, e.g., Q and references therein). 
These include transfer- matrix calculations 0, series expansions ||^ and a scanning Monte 
Carlo method . In two dimensions there is a wealth of exact results for the ordinary surface 
transition {y < y*) from conformal invariance arguments [^^-|T2|,^ and more recently from 
exact Bethe Ansatz calculations More generally, these results have been obtained for 
the 0{n) model, from which the configurational properties of SAWs follow in the n — > 
limit i. 

The situation is not so clear for the special transition. There is a conformal invariance 



result for the thermal scaling dimension X^, which leads to the crossover exponent 



However, the two- variable nature of the special transition poses problems for numerical 
studies, as errors in the estimates of the critical exponents are compounded by errors in the 
location of the critical point {x*, y*). Recent simulations on the square lattice have indicated 
a result significantly larger than the conjectured value 0. 

Here we derive the exact critical couplings {x*, y*) for the polymer adsorption transition 
on the honeycomb lattice, along with the universal critical exponents, from the Bethe Ansatz 
solution of the 0{n) loop model at the special transition. 

Our starting point is the partition function of an 0(n) loop model defined on the 



honeycomb lattice depicted in Fig. 1, 

^ioop = E^V^'^'', (1) 

where the sum is over all configurations of closed and nonintersecting loops. Here P is the 
total number of closed loops of fugacity n in a given configuration. In the limit n this 
reduces to the required SAW generating function, with x the fugacity of a step in the bulk 
and y the fugacity of a step along the surface. Here L is the length of a walk in the bulk 
and Ls is the length of a walk along the surface of the strip. 

The partition function can be conveniently rewritten in terms of the Boltzmann weights 
of the empty vertices. To do this we need to distinguish between three classes of vertices: )— 
and — ( which appear (i) in the bulk and (ii) on the surface, and (iii) ) and ( on the surface. 
For each class we define the weights tf,, ti and t^, respectively. We then consider 

^loop = E^f^'^^'^"'"""^^^^"^^^^''. (2) 

where A/f,, Ni and Ms are the total numbers of vertices (either full or empty) of class (i), 
(ii) and (iii). Apart from harmless normalisation factors, the two partition functions are 
equivalent if U = ti, along with the identification 

x = l/tb, y = l/^Jt^s (3) 

where L = + L^. 

The configurations of the loop model can be mapped to those of a 3-state vertex model 
in the standard way [p^-[T8[|. The allowed arrow configurations and their corresponding 



Boltzmann weights are shown in Fig. 2. Here the phase factors are such that n = s + = 
—2 cos4A. The integrable bulk weights of this honeycomb lattice model are known to follow 
in a particular limit of the Izergin-Korepin model - a more general 3-state model defined on 
the square lattice |]T9|,|l8[. 

For the open boundary conditions of interest here, the integrability of the vertex model 
can be examined in a systematic way by making use of refiection or i^-matrices which satisfy 



the boundary version of the Yang-Baxter equation. In order to do this, we adapted the 



Sklyanin construction of commuting transfer matrices |20| to the present geometry pT|-p3|. 
In particular, we found that the known diagonal reflection matrices for the Izergin-Korepin 
model lead to two integrable sets of boundary weights which preserve the 0(n) symmetry 
2^,|2S|. For each case tb = tj, = 2 cos A. Thus from (0) the critical bulk fugacity is 

1/x* = \/2±V2^, (4) 



which is the well known bulk critical value [|16|. The SAW point occurs at A = vr/S, where 
X* = I/J2 + V2 = 0.541 196 . . .. 



The two integrable sets of boundary weights are ||22 



,., sin2A , cos2A 

(A) ts = —— and (B) = -. 5 

sm A cos A 

Thus from (|]) case (A) gives the critical surface fugacity y* = x* and corresponds to an 

integrable point on the ordinary transition line. However case (B) is new, and corresponds 

to the special transition, with 

y: = {2-n)-'/\ (6) 

in the so-called dilute phase < A < 7r/4 (—2 < n < 2). In contrast the surface coupling 
is complex- valued in the dense phase (vr/4 < A < 7r/2). We thus confine our attention 
here to the dilute region applicable to the adsorption transition. At = we have = 
2^^/^ = 0.840 896 . . . This exact result should prove to be a valuable benchmark for future 
numerical studies of the adsorption transition. For n = 1 the special transition is located at 
y* = 1 which corresponds to infinitely strong surface couplings, as expected. Our result is 
also consistent with a recent argument that although there is no special transition for the 



Ising model, there is a special transition in the geometrical 0{n) model for n > 1 We 
see from that the critical coupling diverges at = 2, and becomes complex for n > 2. 

The central charge c and scaling dimensions Xi defining the critical behaviour of the 
model follow from the dominant finite-size corrections to the transfer matrix eigenvalues 
27|-p9|. The central charge follows from the free energy per site, Jn = log Ag, via 



fN^f^ + J^ + ^,- (7) 

Here fg is the surface free energy and C = 73/2 is a lattice-dependent scale factor. The 
scaling dimensions are related to the inverse correlation lengths via 

= log(Ao/A.) ^ 27rCX,/iV. (8) 

A set of scaling dimensions of interest appear in the so-called watermelon correlator, 
which measures the geometric correlation between ^ nonintersecting SAWs tied together at 
their extremities x and y, which for surface critical phenomena are near the boundary of 
the half-plane [|T2|. It has a critical algebraic decay, 

G,{x-y) ~ \x-y\-^''^. (9) 

These scaling dimensions are associated with the largest eigenvalue in each sector of the 
transfer matrix. In particular, the spin-spin correlation is related to Xi. On the other hand, 
the thermal scaling dimension X^, corresponding to the energy-energy correlation, is related 
to an excitation in the largest sector of the transfer matrix. Given Xi and X^, the surface 
critical exponents follow from 



Xi = 7711/2 = /?i/z/ = [2 + (7 - 27i)/z/]/2, (10) 
X, = 1 - ct>/u. (11) 

In terms of the more standard parameter where 4A + vr^f = 27r, the bulk exponents 7 and 
V are given by 



g 

V 



l={^,g\\lg)v. (12) 



For case (B), at the special transition (a;*, y*), we have obtained the Bethe Ansatz solution 
for the eigenvalues of the transfer matrix 



^ _ -jll sinh(Mj + i 3A/2) sinh(Mj - i 3A/2) 

sinh(Mj + i A/2) sinh(Mj - i A/2) ' ^ ' 



where the Uj follow as roots of the Bethe Ansatz equations 



cosh{uj — i A/2)-|2 rsinh(Mj — i A/2) sinh('Uj — i 3A/2) 



N 



cosh.{uj + i A/2)J ^smh.{uj + i A/2) sinh(Mj + i3A/2) 
-pi sinh(Mj — Ufc + i A) sinh(uj + Uk + i A) sinh(Mj — Uk — i 2A) sinh(uj + Uk — i 2A) 
sinh(Mj — Uk — iX) sinh(Mj + — i A) sinh(Mj — Uk + i 2A) sinh(Mj + Uk + i 2A) 

Here is the width of the strip (e.g., = 8 in Fig. 1) and m labels the sectors of the 
transfer matrix, with m = N for the largest eigenvalue Aq. A more convenient sector label 
is £ = A^ — m. 

The Bethe Ansatz equations differ from those obtained for case (A) in the squared 
pref actor on the left hand side, which have the cosh functions replaced by sinh. This change 
is sufficient to alter the finite-size corrections to the eigenvalues and thus the operator 
content. For case (A), at the ordinary transition, the Bethe Ansatz roots for the largest 
eigenvalue are uniformly distributed along the real (positive) axis |13|. In contrast, at 



the special transition we find that the root distribution includes an elementary 1-string 
excitation, located at mi ~ i (^^ + |^ . To derive the central charge via (]^) we thus adopt the 
analytic method which avoids the explicit manipulation of root densities. Details of the 
calculations will be presented elsewhere. The bulk free energy foo is as derived previously 
17| , p!3[| , while the surface free energy fg differs from that at the ordinary transition ||13[. The 



result is cumbersome and rather unilluminating in the present context. However, at n = 0, it 
reduces to fs = -2 log(l+V2). In fact at this point we observe that Aq = (2+^2)^/(1+^2)^ 
exactly. 

For the central charge, we derive the same result, 

c=l-6{g-l)yg, (15) 



as for the ordinary transition [|1^,|T^,|1^] . We find that the thermal scaling dimension is 
associated with an elementary 2-string excitation, with 

X, = - - 1, (16) 
9 
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in agreement with the conformal invariance result [0. This result has more recently been 
obtained from a thermodynamic Bethe Ansatz calculation from a conjectured boundary 
5- matrix [|26[. From ( p^Tj) and (|1^ the result (^6|) leads to the crossover exponent = ^• 

The root distributions for the eigenvalues defining the the geometric scaling dimensions 
are again real, with 

9-{g-lY 



Xe = lgii+iy-Ui+l) + 



^9 



(17) 



where £=1,2,... These dimensions are to be compared with those at the ordinary transition. 



where the result [12,13 



X, 



\ge + l{g-l% 



(18) 



follows from /i^+i^i in the Kac formula [12 



For the adsorption transition, (|17|) gives 



^9 



(19) 



U{i - 2) 



(20) 



at n = ((7 = I). The first two values are Xi = and X2 = X^ = ^, with the exact 
exponents for the two-dimensional polymer adsorption transition following from ([T0|) and 
(|ll|). In particular, V\\ = ~j2 susceptibility exponent 71 = ||. Guim and Burkhardt 

[0[ originally noted that these were possibly the exact values, as their finite-size scaling 
estimates for Xi and X2 were compatible with X^ = in the Kac formula. Indeed we 

confirm that the more general 0{n) result (p!7| ) agrees with /i£+i,3. Note that the thermal 
dimension (p!6| ) also belongs to the same family of scaling dimensions, since X^ = /ii 3. We 
believe that our results exhaust the complete operator content of the 0{n) model at the 
special transition. 

More recently the result has been conjectured to be correct by Fendley and Saleur 

[|26| who argued that the boundary operator $^+1 3 propagates down the strip at the special 



point. In general, our exact results lend further weight to their claim that the spin degrees of 
freedom of the Kondo problem can be considered as the n = 2 limit of the special transition 
of the 0(?7,) model [p^ . 

So far we have only considered the 0(n) model with boundary conditions that are sym- 
metric with the left and right boundaries of the strip. More generally it is possible to obtain 
the Bethe ansatz solutions of the 0(n) model with non-symmetric boundary conditions 
by using suitable choices of the reflection matrices. In this way we expect to test recent 



conformal invariance results for the 0(n) model with mixed boundary conditions ||3T 
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FIGURES 




FIG. 1. The open honeycomb lattice. 
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FIG. 2. The allowed arrow configurations and corresponding Boltzmann weights for (a) bulk 
and (b) surface vertices. 
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